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ABSTRACT 

Double white dwarf binaries in the Galaxy dominate the gravitational wave sky and would 
be detectable for an instrument such as LISA. Most studies have calculated the expected 
gravitational wave signal under the assumption that the binary white dwarf system can be rep- 
resented by two point masses in orbit. We discuss the accuracy of this approximation for real 
astrophysical systems. For non-relativistic binaries in circular orbit the gravitational wave 
signal can easily be calculated. We show that for these systems the point mass approxima- 
tion is completely justified when the individual stars are axisymmetric irrespective of their 
size. We find that the signal obtained from Smoothed-Particle Hydrodynamics simulations of 
tidally deformed, Roche-lobe filling white dwarfs, including one case when an accretion disc 
is present, is consistent with the point mass approximation. The difference is typically at the 
level of one per cent or less in realistic cases, yielding small errors in the inferred parameters 
of the binaries. 
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1 INTRODUCTION 

At low frequencies (mHz), millions of double white dwarf bina- 
ries in the Galaxy are expected to dominate the gravitational wave 
(GW) sky. At the lowest frequencies they form an unresolved fore- 
ground, while at frequencies above several mHz, thousands of 
sources would be individually detectable for an instrument such 
as the Laser Interferometer Space Antenna, LISA lEvan s~et al.| 
T987l|Lipunov et al.| 1987] |Hils et al.|1990||Nelemans et al.|200i] 
2004 i or eLISA/NGO l |Amaro-Seoane et al.||2012| >. These bina- 
ries come in two flavours: detached systems and semi-detached 
(mass-transferring) systems that are know as AM CVn systems (see 
Solheim 2010 Marsh 2011 for reviews). Several known binaries 
should be detected by LISA within the first weeks of operation 
and are known as verification binaries ( |Stroeer & Vecchio] [2006| 
|Roelofs et al.|2007|[2010l|Brown et al.|2011> . By measuring their 
gravitational wave amplitude and frequency (evolution), the type 
of astrophysical source and its parameters can be determined (e.g. 
|Cutler & FIanagan|1994l|Littenberg|201 l||Blaut|201 1) . In all these 
calculations the gravitational wave signal was determined under the 
assumption that the binary white dwarf system can be represented 
by two point masses in orbit, even for the tidally deformed stars 
in semi-detached binaries. The goal of this study is to determine 
the accuracy of this assumption. In section [2] we will discuss our 
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method of calculating the GW signal. In section [3] we will give 
an algebraic view on the assumption of using point masses and in 
section [4] we will calculate the GW signal from smoothed-particle 
hydrodynamics (SPH) simulations of AM CVn stars. In section [5] 
we discuss the conclusions of this study. 



2 GRAVITATIONAL WAVES FROM ARBITRARY 
SOURCES IN CIRCULAR MOTION 

We calculate the GW wave signal from a collection of point par- 
ticles with arbitrary coordinates and masses, all rotating about a 
fixed point with the same angular speed uj. Since the stars do not 
move at highly relativistic speeds we use linearised general relativ- 
ity. In linearised general relativity the trace reversed metric for any 
non-relativistic, far away source is given by (e.g. Rind ler|2001) 



hi 



-2G f d 



c A R / dt 2 



-pXiXjdV 



-2G 



E 



dt 2 



(1) 

where R is the distance from observer to the source taken as an 
average over the distance to all source points. m a and Xi a are the 
mass and coordinates of points in the source. J? 
the quadrupole moment of the source. In the transverse-traceless 
(TT) gauge for a wave travelling in the z-direction, the wave has 
only two degrees of freedom left. They manifest themselves as so 
called + and x polarisations that can be measured by a detector. 
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The wave metric then becomes (Rindlerj200TJ 



hij 




h+ and h x can be obtained by (Price & Wang|2008| > 



(2) 



(3) 



If we now take N particles and assign them masses m a and posi- 
tions in polar coordinated r a and (f> a = cut + 8 a , we can derive 
a general expression for h+ and h x ■ Also taking into account the 
inclination angle i for the relative orientation of the source, we get: 



c 4 R 



(1 + cos i) (Si cos 2ut - S 2 sin 2wt) (4) 



h x = — ~Tf^ cos i ($2 cos 2«t — Si sin 2ut) , (5) 



where 



Si = m a rj cos 28 a 
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These expressions depend only on the given initial coordinates and 
can easily be calculated numerically. Because h+ and h x have a 
cosine and a sine term with different amplitudes it is useful to define 
the average strain amplitude as: 
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(8) 



As we will later determine the accuracy of the point mass approx- 
imation, it is worth estimating the error we make by using the 
linearised theory for these systems. T he first Post-Newtonian cor 
rection is proportional to (v/c^ 2 



(e.g. 



Misner, Thorne, & Wheeler 



|1973|Blanchet et al.| 19 95 ) which for double white dwarf systems 
with orbital period of several minutes is of order 10 . 



3 PARALLEL AXIS THEOREM FOR GRAVITATIONAL 

WAVES 

Using equations l|4j and {5} we look at what we can say about the 
two point mass approximation algebraically. The parts of the equa- 
tions that depend on the configuration of the system are Si and S2 ■ 
For these expressions something very similar to the parallel axis 
theorem for moment of inertia can be expressed. Consider a body in 
the xy-plane with its centre of mass at the origin. Its mass distribu- 
tion can be described by point masses with coordinates x a ,y a and 
mass m a with respect to its own centre of mass. So Si and S2 are 
obtained from Eqs. |6|7| We call these Si,spmning and S2, spinning, 
because these are the S's that would arise if the body was spin- 
ning around its centre of mass. We now change the coordinates of 
our origin to another point in the plane such that all point masses 
get new coordinates x' a = x a — xcm and y' a — y a — ycu- All 
coordinates denoted with subscript CM refer to the location of the 
body's centre of mass in the new frame and we assume the body is 
in circular motion around the new origin. 

To do the translation we first transform to Cartesian coordi- 
nates, translate and then transform back to polar coordinates. For 



the transform we use 

Xn 



tan(# Q ) for x a ^ 

tan(# Q + 7r) for x a < 



(9) 



The mirror symmetry, cos(2#) = cos(2((9 + n)) and sin(26 l ) = 
sin(2((9 + 7r)) , relieves us from having to split the sum over par- 
ticles into separate sums for positive and negative x a , allowing the 
new S to be written as the old S plus an additional term. The re- 
sulting S's that determine the GW radiation are: 



Si = S 



1 ,spinnin| 



5*2 = Sj.s 



+ ^2 m a r 2 CM cos 26>cm 

a 

+ ^2 m a r 2 CM sin 2#cm, 



(10) 



(11) 



i.e. a first term denoting the original Si. spinning and S2, spinning and 
a second displacement term that sums all mass in one point, so the 
change of coordinates has contributed only the effect of a point 
mass at the body's centre of mass. 

So far we have considered the case in which there is one body 
rotating around an arbitrary point. In the case of a binary system, 
the procedure can be repeated for the second star and the xcm , Vcm 
for each star now refer to the distance of the centres of the two stars 
to the system barycentre. 

We thus can conclude that, when determining h, every body 
with Si, spinning = S2, spinning = can be considered as a point 
mass without implications to the result. These are bodies that do not 
radiate GW if they were only spinning around an axis through their 
own CM. In other words: bodies without a quadrupole moment, 
like axisymmetric spheres or disks. This result even holds when the 
body overlaps with the point it is rotating about, if that is physically 
possible. So only asymmetries in the stars that form a binary GW 
source may lead to deviations from the result as calculated using a 
source consisting only of two point masses. 



4 NUMERICAL CALCULATION OF GRAVITATIONAL 
WAVES FROM SPH SIMULATIONS 

For semi-detached binaries, we know they do not consist of spheri- 
cal stars. The contribution of the accretion disk will depend on if it 
is circular or not and on how the mass is distributed over the disk. 
The contribution of the donor will also depend on the mass distri- 
bution over its shape. We therefore need to look at the GW signals 
of more realistic mass distributions. 



4.1 The simulations 

Roche lobe filling stars and accretion disks are not completely sym- 
metric around their centres of mass. Using equations (|4| and |5j we 
can calculate the contribution this has to the gravitational wave sig- 
nal if we know their mass distributions. For this study we use a set 
of SPH simulations of double white dwarfs at the onset of mass 
transfer (taken from |Dan et al. 201 1 2012 1, hereafter referred to as 
RocheSPH. We sample the different mass combinations as shown 
in Table [T] An example of one star of an SPH simulation is shown 
in Fig.[T|in such a way that the non-axisymmetric SPH particles are 
shown more prominently. In addition we have used an SPH simula- 
tion of an accretion disc in an ultra-compact binary, kindly provided 
to us by Prof. Matt Wood (based on Woo d, Thomas, & Simpson] 
|2009| >. The latter is a Mdonor/A/accretor = q = 1/10 system with 
an accretion disk around the accretor formed by adding mass at 
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Figure 1. Projection of SPH particles onto the orbital plane of one of 
the stars in a double white dwarf system with Mj onor = iWaccretor = 
0.2Mq ,q = l. The colour coding accentuates the asymmetry in the mass 
distribution that is responsible for the deviation of the gravitational wave 
signal from that of two point masses. For each ring around the center of 
mass of the star the ratio of the mass on the left (x coordinates smaller 
than xcm) t° that on the right(x coordinates larger than xcm) is calcu- 
lated. Light grey represents values of unity (smoothed for clarity) and black 
values below 0.5. 
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Figure 3. Ratio of the GW strain amplitude (equation [8j from the 
RocheSPH systems (Fig.[TJ and corresponding the two point mass systems 
as a function of the mass ratio, for different donor masses. See also Table[T] 



Table 1. Ratio of GW strain amplitude of SPH systems and two point 
mass systems. The mass of the accretion disk in DiskSPH was taken to 

be 10 ^-^^accretor ■ 
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Figure 2. Projection of SPH particles onto the orbital plane of a double 
white dwarf system with q = 1/10, the donor and accretor are represented 
as point masses (the bigger dots) and an accretion disk is present around 
the accretor. The coordinates are in units of the separation a between the 
two stars. In this SPH simulation particles were continuously added at the 
inner Lagrange point of the Roche potential. These particles have formed 
the accretion disk around the accretor. This is a picture of the 390th orbit 
after the mass transfer started. 
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the inner Lagrange point of the Roche potential with the two stars 
treated as point masses (Fig. [2J, from here onwards referred to as 
the DiskSPH. We assume the mass in the disk is 10 _J of the mass 
of the accretor. 

All simulation results are snapshots at some point in the evo- 
lution of the binary. Because of this, we can not take all movement 
of the system into account, so we assumed that these mass distribu- 
tions rotate with fixed angular speeds around their centres of mass. 
Now equations [4] [5] and [8] can be used to calculate the GW strain 
amplitude. 

4.2 The results 

The differences of GW strain amplitude for the DiskSPH and all 
RocheSPH systems compared to two point systems are listed in Ta- 
ble [I] The DiskSPH results shows a deviation at the 10 -5 level, 



while the RocheSPH results, depending on the mass of the de- 
formed donor star and the mass ratio q, range between 0.2 and 1.3 
per cent. The results of the RocheSPH calculations as function of 
mass ratio q and donor mass are shown graphically in Fig. [3] The 
deviations are largest for equal mass systems, but more interesting, 
our results clearly show the fact that lower donor masses are more 
deformed than more massive donors at the same mass ratio. 

To explore the influence of an accretion disk in somewhat 
more detail, we varied the mass of the accretion disk compared 
to the mass in the two stars (Fig. H). As can be expected, the de- 
viation from the point mass approximation scales linearly with the 
disk mass, which has to be unrealistically high (1 per cent of the 
total mass) in order to get near the effect of the deformed donor 
star discussed above. The deviation for a system with Roche-lobe 
filling donor plus an accretion disk is approximately the sum of the 
RocheSPH and DiskSPH result (see above) and thus for realistic 
disk masses is dominated by the deformation of the donor. 
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Figure 4. Ratio of the GW strain amplitude (equation[8| from the DiskSPH 
system (Fig. [5) and the two point mass system. The x-axis shows ratio of 
disk mass to the accretor plus disk mass. The donor is considered to be a 
point mass, so at m<ji s k = the system is equal to the two point system. 



5 CONCLUSIONS 

We have calculated the deviation of the gravitational wave signal of 
finite size and non-spherical white dwarfs to that of point masses, 
which is usually assumed. For any finite size, non-relativistic ax- 
isymmetric body in circular orbit, the result is exactly the same. 
For semi-detached white dwarf binaries, we find that an accretion 
disk of reasonable mass changes the gravitational wave signal at 
the level of to 10~ 3 , small but still significantly larger than 
errors due to the neglect of post-Newtonian corrections in the cal- 
culation of the signal. Deformations due to filling the Roche lobe 
of semi-detached binaries increase for mass ratios closer to unity 
for fixed donor mass and are distinctly stronger at fixed mass ra- 
tio for lower mass donors and can in the most extreme cases be of 
order 1 per cent. This in principle will change the frequency evolu- 
tion and the accuracy with which the parameters can be determined 
(e.g. |Blanchet et al.|1995)|Cutler & Flanagan|1994| >, so calculations 
in which accuracies of better than one per cent are needed should 
take the finite size into account. Also, the different strength of the 
gravitational wave angular momentum losses will affect the mass 
transfer rate and stability of the mass transfer in semi-detached sys- 
tems (e.g. Marsh et aL 2004). However, the level of deviation we 
found shows that the calculations presented in the literature on the 
expected signals of (verification) binaries for LISA, are essentially 
unaffected: for monochromatic sources, the amplitude of the sig- 
nal will be slightly different, but indistinguishable from sources 
with slightly higher (chirp) masses and/or smaller distances, prop- 
erties that are much more uncertain than one per cent. For systems 
with measurable period derivatives, the degeneracy between (chirp) 
mass and distance is broken, but the deviations we found can still 
only be detected if there is independent extremely accurate mea- 
surement of the (chirp) mass of the binary. For a LISA-like de- 
tector such independent mass estimates, if available, are typically 
accurate at the 10 per cent level at best (e.g. Littenberg 201 1 ). The 
detailed evolution (and possible merger of the system) will be af- 
fected as well, but the calculations of this phase (e.g. |Dan et aT] 
|2012[ l already take the finite size of the stars into account. We there- 
fore conclude that in the majority of cases, the use of the point mass 
approximation is well justified. 
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